REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and 
Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person 
shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR 
FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

July  2015  Briefing  Charts 


4.  TITLE  AND  SUBTITLE 

Theoretical  Studies  of  Nanoclusters  (Briefing  Charts) 


3.  DATES  COVERED  (From  -  To) 
July  20 15- July  2015 


5a.  CONTRACT  NUMBER 

In-House 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/RQRP 
10  E.  Saturn  Blvd. 

Edwards  AFB,  CA93524-7680 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/RQR 
5  Pollux  Drive 

Edwards  AFB  CA  93524-7048 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited. 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 

Q188 


8.  PERFORMING  ORGANIZATION 
REPORT  NO. 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

AFRL-RQ-ED- V  G-2 0 1 5-2 54 


13.  SUPPLEMENTARY  NOTES 

Briefing  Charts  presented  at  Air  Force  High  Performance  Computing  User  Forum;  Dayton,  OH;  23  July  2015.  PA#15341. 


14.  ABSTRACT 

Briefing  Charts 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

Unclassified 

Unclassified 

Unclassified 

17.  LIMITATION 
OF  ABSTRACT 


18.  NUMBER 
OF  PAGES 


19a.  NAME  OF 
RESPONSIBLE  PERSON 

Jerry  Boatz 


19b.  TELEPHONE  NO 

(include  area  code) 

661-275-5364 


Standard  Form 
298  (Rev.  8-98) 

Prescribed  by  ANSI 


Theoretical  studies  of  nanoclusters 

23  Jul  2015 


Principal  Research  Chemist 
Aerospace  Systems  Directorate,  RQRP 
Air  Force  Research  Laboratory 

This  briefing  contains  information  up  to: 


Distribution  A:  Approved  for  public  release;  distribution 
unlimited.  PA  Clearance  Number  15341 


1.  Introduction 

-  Modeling  and  simulation  methods 

-  objectives,  payoffs,  approach 

2.  Aluminum  nanoparticles 

-  energetic  additives  for  propellants,  explosives 

-  “hybrid”  fuels  consisting  of  nanoparticles  suspended  or 

dissolved  in  ionic  liquids 

-  surface-functionalized  NPs  can  be  dispersed  in  variety  of  liquids 

hydrocarbons,  ILs,  polar  solvents,  etc. 

3.  Core-shell  nanoclusters 

-  energetic  additives  for  propellants,  explosives 

-  gas  generators 

-  biocidal  defeat  agents 

4.  Summary  and  conclusions 
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Objectives 


*  Support  and  streamline  the  discovery,  synthesis, 
and  characterization  of  new  energetic  materials  for 
chemical  propulsion,  explosives,  gas  generators, 
etc. 

-  predictions  that  are  accurate,  relevant,  and  timely. 

*  Utilize  modeling  and  simulation  as  computational 
tool  to 

-  identify  suitable  target  compounds 

-  “what-if”  scenarios,  trial-and-error  explored  via  M&S 

-  explore  possible  synthesis  routes 

-  confirm  successful  synthesis 

-  provide  detailed  analysis  of  chemistry  &  mechanisms 

*  Expand  scope  and  complexity  of  problems  which 
can  be  addressed  via  M&S  &  improve  M&S  reliability 
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Payoffs 


*  Focus  experimental  efforts  in  the  most  promising 
directions 

*  Identify  technical  “dead  ends”  early  in  the  propellant 
development  cycle 

*  Provide  fundamental  insight  into  observed  behavior 
of  materials 

*  Enable  rational  design  approaches  in  development 
of  advanced  materials 


Distribution  A:  Approved  for  public  release;  distribution  unlimited. 


4 


Technical  Approach 


•  Utilize  robust,  physics-based  M&S 

•  Quantum  chemical  methods  (molecular  quantum  mechanics) 

•  Predictive  capability 

•  no  parameterization,  empirical  fitting,  etc. 


•  Leverage  fundamental  R&D 

•  AFOSR  support  for  computational  chemistry 

-  e.g.,  development  of  “fragment”  methods  for  large-scale 
computations 

•  NRC  postdocs,  senior  fellowships,  summer  faculty  fellowships 


•  Leverage  DoD  HPC  resources 

•  -100M  cpu  hours  /  year 

•  HPC  frontier,  HPC  HASI,  and  PETTT  pre-planned  projects 

•  Institutes  (e.g.,  Multi-Scale  Reactive  Modeling) 
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Aluminum  nanoparticles 


•Al  nanoparticles  (NPs)  are  of  interest  as  energetic  ingredients  in 
explosives  and  propellant  formulations,  due  to  high  energy 
density,  enhanced  burn  rates,  etc. 


•  Efficient  production  of  Al  nanoparticles  via  ball  milling  is  obtained 

using  NH3,  CH3NH2,  orCH3CN  as  milling  agents. 

•  Milling  agents  decompose  on  NP  surface  to  produce  gaseous 

products  and  surface-bound  species. 

•  Milling  agents  can  also  passivate  NP  surface  against  oxide  layer 

formation 
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Ball  milling  gaseous  byproducts 


Milling  agent/ 
Gaseous  products 

NH3 

ch3nh2 

ch3cn 

H2  1 

V 

V 

V  I 

CH3NHCH3  y 


ch2nh 

V 

ch4  V 


ch3ch3 

V 

Can  theory  explain,  for  example  in  the  case  of  NH3  as  the 
milling  agent,  the  predominant  formation  of  H2  as  well  as  the 
absence  of  other  stable  NxHy  species  (N2,  N2H2,  N2H4?) 
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Quantum  chemical  calculations 


•  DFT  calculations  (M06/6-311++G(d,p))  used  to  identify  surface 

species  and  reactions  leading  to  their  formation. 

•  Al80  nanocluster  used  as  NP  model. 

-  contains  bulk-like  “core”  partially  surrounded  by  surface  layer 

-  surface  is  atomically  “rough”  -  realistic  representation  of  actual  NP 
surface? 

•  Reaction  enthalpies  and  barriers  for 

-  chemisorption  of  NH3 

-  fragmentation  of  chemisorbed  NH3 

-  formation  of  H2 

•  direct  elimination  of  H2  from  chemisorbed  NH3 

•  dissociative  recombination  of  chemisorbed  H  atoms 

•  dissociative  elimination  of  H2  from  adjacent  chemisorbed  NHX  (x=1-3) 

-  formation  of  N2,  N2H2,  N2H4  (not  observed  in  experiments) 
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NH3  chemisorbed  on  Al80 
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H3N:Alao 


NH2-AIs0-H 


AHrxn  =  -2  kcal/mol 
AH*  =  +31  kcal/mol 


•• 


+19.6  (+16.7) 
vimag=1549i  cm-1 
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Fragmentation  of  NH3  on  Al80 
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Formation  of  H2 
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Reaction  coordinate  (bohr-amu1/2) 
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a)  ^  Formation  of  N2,  N2H2,  N2H4 


Core-shell  nanocluster  synthesis 

Core-shell  nanoclusters  such  as  SiAln,  NinAlm,  Aln(CuO)mJ  etc.  may  be  useful 
ingredients  in  propellants  and  explosives 

-  higher  energy  densities  than  organics  (~  3x  RDX) 

_ -  some  are  resistant  to  surface  oxidation  (i.e.,  “magic  clusters”) 


Helium  droplet 
experiments  at 
AFRL/RW 


Can  core-shell  nanoclusters  be  formed  under  cryogenic  conditions  (i.e.,  in 
helium  droplet  experiments)  via  stepwise  condensation;  i.e.,  what  are  the 
energy  barriers  (if  any)  to  stepwise  addition  of  atomic  Al? 

_ SiAln  +  Al  -►  SiAln+1 .  [Aln]-  +  Al  -►  [Al]~ _ 
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CuxMgy  core-shell  nanocluster 

inversion 


In  helium  droplet  experiments,  Mg  atoms  were  captured  in  first  pickup 
cell,  followed  by  capture  of  Cu  atoms  to  form  CuxMgy  core-shell 
nanoclusters.  However,  scanning  transmission  electron  microscopy 
(STEM)  measures  show  cluster  inversion  occurred  to  produce  MgyCux(!) 


a)  copper  atoms 

b)  magnesium  atoms 

c)  oxygen  atoms 

d)  composite  image 


Distribution  A:  Approved  for  public  release;  distribution  unlimited. 


15 


1 .  Structure  of  Mg30  cluster 
was  fully  optimized. 

2.  Two  Cu  atoms  were 
placed  on  opposite  sides  of 
Mg30  and  structure 
reoptimized. 

3.  Distance  between  Cu 
atoms  was  decreased  in 
steps  of  0.25  A,  held  fixed, 
and  remaining  DOF 
reoptimized. 

4.  Total  energy  plotted  as 
function  of  fixed  Cu-Cu 
distance. 


-16  i  ■  ■  i  ■  i  ■  i  ■  i 

5  5.5  6  6.5  7  7.5  8  8.5  9  9.5  10 

r  (Cu-Cu)  (A) 


DFT  calculations:  B3PW91/aug-cc-pcoCVTZ(-PP)  level 
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Cu2Mg30  vs.  Mg30 
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Need  for  HPC 


•  Individual  job  characteristics 


•  Al80  +  NH3  (per  individual  system) 

•  Structure  optimization:  600  -  1200  cores  x  200  hours  wall  time 

=  120-240K  core-hours 

•  Vibrational  frequencies:  2500  -  5000  cores  x  200  hours  wall  time 

=  500K-1M  core-hours 

•  Intrinsic  reaction  coordinate:  600-1200  cores  x  400  hours  wall  time 

=  240-480K  core-hours 

•  Cu2Mg30 

•  Structure  optimization:  2400-3200  cores  x  100  hours  wall  time 

=  240-300K  core-hours 

•  Vibrational  frequencies:  2400-3200  cores  x  300  hours  wall  time 

=  720-960K  core-hours 
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Summary  and  Conclusions 


•  Production  of  Al  nanoparticles  via  ball  milling 

•  NH3-assisted  ball  milling  of  Al  powder  efficiently  produces  Al  NPs,  with  H2  as 
the  predominant  byproduct. 

•  Surface  reactions  of  NH3  on  Al80  have  been  modeled  using  DFT 

•  NH3  chemisorbs  to  Al80  with  binding  energies  of  8-1 6  kcal/mol 

•  Surface  fragmentation  of  NH3  to  form  chemisorbed  NH2  +  H  is  slightly  exothermic 
but  has  a  barrier  of  30  kcal/mol. 

•  Elimination  of  H2  from  adjacent  chemisorbed  H  and  NH3  is  slightly  endothermic  but 
has  a  barrier  of  only  12  kcal/mol. 

•  Formation  of  N2,  N2H2,  N2H4  from  recombination  of  chemisorbed  N,  NH,  and  NH2 
respectively,  is  endothermic  by  ~42,  103,  70  kcal/mol. 

•  Calculations  are  consistent  with  observation  of  H2  as  the  predominant 
byproduct  and  minimal  amount  of  N2. 


•  Mg/Cu  core-shell  nanoclusters 

•  Helium  droplet  experiments  show  inversion  of  CuxMgy  clusters  to  MgyCux. 

•  Cu  atoms  diffusing  into  Mg30  and  vice-versa,  have  been  modeled  using  DFT. 

•  Estimated  barrier  for  Cu  atoms  to  migrate  into  Mgn  is  <  1  kcal/mol. 

•  Estimated  barrier  for  Mg  atoms  to  migrate  into  Cun  is  6  kcal/mol. 

•  Calculations  are  consistent  with  observed  Cu/Mg  inversion. 
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Planetary  ball  milling 
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Technical  Progress:  EILs 


•  Can  metal  nanoclusters  (e.g.,  aluminum  or  boron)  form  stable  colloidal  suspensions  in  ionic  liquids  (ILs)  to 
produce  air-stable,  oxidation-resistant,  hybrid  fuels? 

•  Why  do  boron  nanoparticles  (BNPs)  form  stable  suspensions  in  some  solvents  but  not  others? 

■  What  types  of  chemical  interactions  occur  between  solvent/IL  and  BNPs? 


Boron  nanopowder  milled  in 
IL/ethanol  mixture: 

STABLE 

Boron 

nanoclusters 
remain  in 
solution;  form 
a  stable 
colloidal 
suspension 


■  Utilize  M&S  to 

•  understand  the  difference  in  behavior  of  BNPs  milled  in  conventional  solvents  vs.  ILs 
-  optimize  the  long-term  stability  of  IL/BNP  colloidal  suspensions 


Customers: 

RQRP  /  Univ.  of  Utah  /  Univ.  of  Alabama 

Involvement: 

experimental  synthesis  &  data 

Sponsor: 

AFOSR 
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IL/metal  hybrid  fuels 


•  Addition  of  1 0%  aluminum  or  boron  into  each  bipropellant  yields  small 
changes  in  specific  impulse  for  each  fuel  (MMH,  RP1  &  IL) 

•  Aluminum  gives  slightly  more  positive  improvement  than  boron 


•  Aluminum  &  Boron  nano-particulate  significantly  increases  volumetric 
impulse  (energy  density)  for  all  3  fuels 

•  MMH/NTO  does  maintain  Isp  advantage  over  both  hydrocarbon  (RP1)  and 
dicyanamide-based  IL  fuels  -  due  to  lower  MW  exhaust  species 

•  Volumetric  impulse  of  IL  fuel  without  metal  hybridization  is  good 

•  Same  as  that  of  boron-bearing  MMH  and  hydrocarbon  fuels 

•  Greater  than  aluminum-bearing  MMH  and  hydrocarbon  fuels 

•  IL  fuel  provides  highest  volumetric  impulse  of  all  three  fuels  when 
hybridized  with  metal 

•  Boron  &  aluminum  appear  to  provide  similar  predicted  IL  fuel  performance  (metal 
held  at  1 0  wt%  of  totai  propellant) 

•  Aluminum  is  more  effective  (i.e.,  lower  concentration  of  metal  required  in  fuel 
component)  compared  to  boron  to  achieve  same  performance 

•  For  IL  fuel,  21 .7  wt%  Al  for  14.2  lbf-s/in3,  versus  27.0  wt%  B  for  same 
performance 
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IL/metal  hybrid  fuels 


Computational  model  -  why  B80?? 


-  Calculations  predict  most  stable  forms  of  Bn  for  n  <  20  are  quasi-planar 

-  Most  stable  form  of  B20  is  a  ring 


Table  I.  The  MP2/6-311G*  optimized  B20  geometries  and  computed  relative  energies  (in  eV)  of  the  eight  B20  isomers  at  different  theoretical  levels. 


Symmetry 

Qv 

C2v 

c2 

Q 

c, 

cs 

c. 

S4 

CCSD(T) 

0 

0.72 

1.46 

1.87 

1.97 

2.31 

2.80 

3.45 

MP2 

0 

1.13 

1.07 

0.94 

1.89 

1.47 

1.16 

3.32 

PBE 

0 

0.67 

1.51 

1.69 

2.26 

2.52 

2.60 

3.80 

TPSS 

0 

0.68 

1.66 

1.96 

2.17 

2.46 

3.00 

3.35 

TPSSh 

0 

0.79 

1.65 

1.93 

2.12 

2.44 

3.09 

3.23 

PBEO 

0 

0.96 

1.50 

1.68 

2.13 

2.47 

2.90 

3.53 

mPWlPW91 

0 

0.99 

1.75 

2.15 

2.37 

2.74 

3.74 

3.75 

M06-2X 

0 

1.11 

1.03 

0.93 

1.61 

1.98 

1.13 

2.71 

B3LYP 

0 

0.99 

3.25 

4.24 

3.82 

4.39 

4.80 

5.34 

BLYP 

0 

0.75 

3.38 

4.39 

4.04 

4.58 

4.78 

5.63 

Fengyu  Li,  Peng  Jin,  De-en  Jiang,  Lu  Wang,  Shengbai  B.  Zhang,  Jijun  Zhao,  and  Zhongfang  Chen,  J.  Chem.  Phys.  136,  074302  (2012) 
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IL/metal  hybrid  fuels 


Computational  model  -  why  B80?? 


-  IL  interactions  with  small  3D  clusters  (e.g.,  B12)  inconsistent  with  larger  clusters 


B12  +  [MAT][DCA] 

B80  +  [MAT][DCA] 

Eint  =  104.5  kcal/mol 

Eint  =  51.5  kcal/mol 

Distinct  differences  in 

-  interaction  energies  (2x) 

-  binding  of  cation 

-  binding  of  anion 


Distribution  A:  Approved  for  public  release;  distribution  unlimited. 
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Mgn  benchmark  calculations 


Mgn  clusters 

“Closed  shell”  atomic  configuration  [(1s)2(2s)2(2p)6(3s)2]  suggests  that 
weak  dispersion  interactions  will  be  important.  Need  to  consider 

-  core-core  and  core-valence  correlation 

-  correlation  method  (MP2,  CC,  DFT) 

-  “active”  electrons  to  be  correlated  in  MP2,  CC 

-  suitable  DFT  functional  for  larger  Mgn 

clusters  (up  to  n  ~  100) 


Method 

cc-pwCVDZ 

cc-pwCVTZ 

cc-pwCVQZ 

MP2 

23.1 

/ 

3.042 

28.5 

/ 

3.013 

29.1 

/ 

3.011 

CCSD(T) 

16.7 

/ 

3.100 

tbd 

/ 

3.064 

tbd/ 

3.065 

DFT/B3PW91 

26.5 

/ 

3.092 

26.5 

/ 

3.091 

26.4 

/ 

3.092 

DFT/PBE 

34.6 

/ 

3.070 

34.5 

/ 

3.070 

34.3 

/ 

3.070 

DFT/PBEO 

31.5 

/ 

3.078 

31.5 

/ 

3.078 

31.4 

/ 

3.078 

DFT/M06 

30.8 

/ 

3.028 

30.4 

/ 

3.025 

TBD 

/ 

TBD 

D  FT/M  11 

19.2 

/ 

3.134 

TBD 

/ 

TBD 

TBD 

/ 

TBD 

Calculated  binding  energies  used  to  determine  size  of  helium  droplet  needed  for  evaporative  cooling 
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